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ffl generating and ranking initial conformations for each ligand in the set of ligands at the binding 
region using docking techniques ; 



optimizing the preferred binding conformations using annealing molecular dynamics, the 



calculating a binding energy for each ligand in efthe set of ligands in the corresponding 
optimized preferred binding conformations; and 

selecting for each ligand in eaoh-of the set of ligands the lowest calculated binding energy 
in the optimized preferred binding conformations, and outputting the selected calculated binding 
energies as the predicted binding energies for the predicted binding conformations' for each of the 



set of ligands. 

2. (Original) The method of claim I, wherein: 

the binding region is a known binding region defined by the structural information . 

3. (Currently amended) The method of claim 1, wherein: 

« 

the binding region is an unknown binding region; and 

using the structural information for the protein to identify a binding region of the 
ligand in the p rotein comprises predicting a probable binding region based at least in part on the 
structural information. 

4. (Currently amended) The method of claim 3, wherein using the structural 
information to identify a binding region of the protein predicting a probable binding region 
comprises: 

mapping the empty volumes available for ligand binding in the protein to identify one or 
more potential binding regions; 

generating initial conformations for one or more ligands known to bind the protein using 
docking techniques in each of the one or more potential binding regions and scoring a 
preliminary energy function for at least some of the initial conformations; 
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selecting from the initial conformations for each of the known ligands a plurality of best 

* 

conformations in each of the potential binding regions based at least in part on the preliminary 
energy scores; 

» 

and scoring an energy function for each of the best conformations ; and 

identifying the probable binding site based on a spatial location of the test conformations 

5. (Currently amended) The method of claim 4, further comprising: 

before s coring the energy function for e ach of the b e st conformations, optimizing the 
selected best conformations to obtain a set of energy-minimized conformations for each of the 
known ligands in each of the potential binding regions and scoring a second preliminary energy 
function for each of the best conformations ; 

wherein identifying the probable binding site is based on a spatial location of the best 
conformations having the lowest second preliminary energy scores, wherein the energy function 
is scored for each of the energy minimized conformations, 

* 

6. (Currently amended) The method of claim 4, further comprising: 
before scoring the preliminary energy function for each of the best conformations, 

calculating for each of the best conformations a percentage of the ligand surface area buried 
within the protein for the conformation; 

wherein the preliminary energy function is scored only for a subset of the best 
conformations, wherein each of the best conformations in the subset has having a calculated 
percentage of the ligand surface area buried within the protein which exceedsmg-a predetermined 
surface area threshold. 

7. (Cancelled) 

8. (Currently amended) The method of claim 71, further comprisi ng wherein : 
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identifying a plurality of preferred binding conformations includes after s electing the best 
conformations, optimizing the selected best conformations to obtain a set of energy-minimized 
conformations for each of the ligands; 

wherein the preferred binding conformations comprise the energy-minimized 
conformations. 

9. (Original) The method of claim 1, wherein: 

the annealing molecular dynamics includes a full atom force field. 

10. (Original) The method of claim 1, wherein: 

the solvation effects include a continuum description of solvation. 

1 1 . (Original) The method of claim 1 , wherein: 

the solvation effects include a surface-area based solvation model. 

12. (Currently amended) The method of claim 1, wherein: 

calculating a binding energy for each Ugandjn of the set of ligands includes taking the 
difference in the ligand energy in the receptorprotein and in solution. 

9 

13. (Original) The method of claim 1, wherein: 

the binding energy is calculated for a ligand according to a scoring function comprising 
subtracting the free energy of the ligand in water from the energy of the ligand in the protein. 

14. (Original) The method of claim 1, wherein: 

the binding energy is calculated for a ligand according to a scoring function comprising 
subtracting the free energy of the protein and the free energy of the ligand from the free energy 
of the ligand in the protein. 

15. (Original) The method of claim 1 , further comprising: 



f 
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identifying from the set of iigands one or more ligands predicted to have high binding 
affinity based on the calculated binding energy of the ligands in the binding site. 

16. (Original) The method of claim 1, wherein: 

the protein is a globular protein or a transmembrane protein. 

1 7. (Withdrawn) A computer-implemented method for predicting the structure of a 
protein binding site for a protein having an unknown binding site, the method comprising: 

providing structural information describing the structure of a protein having an unknown 
binding site and a set of one or more ligands known to bind to the protein; 

using the structural information for the protein to identify a plurality of potential binding 
regions of the protein; 

generating initial conformations for one or more of the ligands using docking techniques 
in each of the potential binding regions; 

selecting from the initial conformations for each of the ligands a plurality of best 

* 

conformations in each of the potential binding regions and scoring an energy function for each of 
the best conformations; 

identifying the probable binding site based on a spatial location of the conformations 
having the lowest energy scores; and 

outputting structure information describing the three-dimensional structure of the 
probable binding site. 

1 8. (Withdrawn) The method of claim 1 7, further comprising: 

before scoring the energy function for each of the best conformations, optimizing the 
selected best conformations to obtain a set of energy-minimized conformations for each of the 
ligands in each of the potential binding regions; 

wherein the energy function is scored for each of the energy-minimized conformations. 

19. (Withdrawn) The method of claim 17, further comprising: 
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before scoring the energy function for each of the best conformations, calculating for 
each of the best conformations a percentage of the ligand surface area buried within the protein 
for the conformation; 

> 

wherein the energy function is scored only for a subset of the best conformations having 
a calculated percentage of the ligand surface area buried within the protein exceeding a 
predetermined surface area threshold. 

20. (Withdrawn) A computer-implemented virtual screening method for screening a 
ligand library, the method comprising: 

receiving protein structural information describing the structure of a protein; 
receiving ligand structural information describing the structure of a plurality of ligands in 
a ligand library; 

receiving an input specifying a desired number of candidate ligands to be identified in the 
ligand library; 

using the structural information for the protein to identify a binding region of the protein; 
generating a set of initial binding conformations for each of the ligands in the binding 

region; 

calculating an energy function for each of the initial binding conformations and selecting 
for each of the ligands a plurality of the initial binding conformations having the lowest 
calculated energy as a set of best conformations; 

optimizing the best conformations; 

calculating a binding energy for each of the ligands in the corresponding optimized best 
conformations; and 

selecting from the plurality of ligands a set of the desired number of candidate ligands 
having the lowest calculated binding energy in the optimized best binding conformations, and 
outputting the selected set of candidate ligands. 

♦ 

21. (Withdrawn) The method of claim 20, wherein: 
the plurality of ligands comprises at least SOO ligands* 
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22. (Withdrawn) The method of claim 20, wherein: 
the plurality of ligands comprises at least 1 ,000 ligands. 

23. (Withdrawn) The method of claim 20, wherein: 

* 

the plurality of ligands comprises at least 5,000 ligands. 

24. (Withdrawn) The method of claim 20, wherein: 
the plurality of ligands comprises at least 1 0,000 ligands. 

25. (Withdrawn) The method of claim 20, wherein: 
the plurality of ligands comprises at least 50,000 ligands. 

26. (Withdrawn) The method of claim 20, wherein: 
the plurality of ligands comprises at least 100,000 ligands. 

t 

27. (Withdrawn) The method of claim 20, wherein: 

calculating a binding energy for each of the set of ligands includes taking the difference 
in the ligand energy in the receptor and in solution. 

28. (Withdrawn) The method of claim 20, wherein: 

the binding energy is calculated for a ligand according to a scoring function comprising 
subtracting the free energy of the ligand in water from the energy of the ligand in the protein. 

29. (Currently amended) A computational model of a ligand-protein complex for a 

■ 

protein having an unknown binding site, the model comprising: 

a computer-readable memory storing data describing an optimized preferred binding 
conformation for the protein and a ligand known to bind to the protein, the optimized binding 
conformation being generated according to the method of claim 1. 



i 
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30. (Withdrawn) A computational model of a predicted structure for a protein 
binding site for a protein having an unknown binding site, the model comprising: 

a computer-readable memory storing data describing the three-dimensional structure of 

the probable binding site for the protein generated according to the method claim 1 5. 

» 

9 

3 1 . (Currently amended) A computer program product on a computer-readable 
medium for modeling ligand-protein binding interactions, the computer program product 
comprising instructions operable to cause a programmable processor to: 

provide structural information describing the structure of a protein and each ligand in a 
set of one or more ligands; 

use the structural information for the protein to identify a binding region of the protein; 

identify a plurality of preferred binding conformations for each ligand in ef-the set of 
ligands in the binding regio n, the preferred binding conformations being determined by 
generating and ranking initial conformations for each ligand in the set of ligands at the binding 
region using docking techniques ; 

optimize the preferred binding conformations using annealing molecular dynamics, the 
annealing molecular dynamics including solvation effects; 

calculate a binding energy for each ligand in ef-the set of ligands in the corresponding 
optimized preferred binding conformations; and 

select for each ligand in of-the set of ligands the lowest calculated binding energy in the 
optimized preferred binding conformations, and output the selected calculated binding energies 
as the predicted binding energies for the predicted binding conformations for each of the set of 
ligands, 

32. (Withdrawn) A computer program product on a computer-readable medium for 
predicting the structure of a protein binding site for a protein having an unknown binding site, 
the computer program product comprising instructions operable to cause a programmable 
processor to: 
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* 

provide structural information describing the structure of a protein having an unknown 
binding site and a set of one or more ligands known to bind to the protein; 

use the structural information for the protein to identify a plurality of potential binding 
regions of the protein; 

generate initial conformations for one or more of the ligands using docking techniques in 
each of the potential binding regions; 

select from the initial conformations for each of the ligands a plurality of best 
conformations in each of the potential binding regions and score an energy function for each of 
the best conformations; 

identify the probable binding site based on at least one of a percentage surface area of the 
ligand buried in the protein or a spatial location of the conformation having the lowest energy 
score; and 

output structure information describing the three-dimensional structure of the probable 
binding site. 

33. (Withdrawn) A computer program product on a computer-readable medium for 
screening a ligand library, the computer program product comprising instructions operable to 

« 

cause a programmable processor to: 

receive protein structural information describing the structure of a protein; 

receive ligand structural information describing the structure of a plurality of ligands in a 
ligand library; 

receive an input specifying a desired number of candidate ligands to be identified in the 
ligand library; 

use the structural information for the protein to identify a binding region of the protein; 
generate a set of initial binding conformations for each of the ligands in the binding 

region; 

calculate an energy function for each of the initial binding conformations and selecting 
for each of the ligands a plurality of the initial binding conformations having the lowest 
calculated energy as a set of best conformations; 
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optimize the best conformations; 

calculate a binding energy that includes solvation for each of the ligands in the 
corresponding optimized best conformations; and 

select from the plurality of ligands a set of the desired number of candidate ligands 
having the lowest calculated binding energy in the optimized best binding conformations, and 
output the selected set of candidate ligands. 

34. (Withdrawn) A computer-implemented method of generating a pharmacophore, 
comprising: 

providing structural information describing the structure of a protein and a set of one or 

■ 

more ligands known to bind to the protein; 

using the structural information for the protein to identify a binding region of the protein; 

identifying a plurality of preferred binding conformations for each of the set of ligands in 
the binding region; 

optimizing the preferred binding conformations using annealing molecular dynamics, the 
annealing molecular dynamics including solvation effects; 

calculating a binding energy for each of the set of ligands in the corresponding optimized 
preferred binding conformations; 

selecting for each of the set of ligands the optimized preferred binding conformation 
having the lowest calculated binding energy; 

generating a pharmacophore model based at least in part on the selected optimized 
preferred binding conformations, the pharmacophore model defining a pattern of ligand features 
predicted to be required for binding to the protein; and 

outputting data representing the pharmacophore model for use in drug design. 

35. (Withdrawn) The method of claim 34, further comprising: 

using the pharmacophore model as a template to search a chemical information database 
to identify one or more molecules predicted to bind to the protein. 



Applicant : Wely B. Floriano, Nagarajan Vaidehi, Attorney's Docket No. : 066 1 8-607002 / CIT 3 1 92-C 

William A. Goddard, III 

Serial N . : 10/010,725 

Filed : November 30, 2001 

Page : 11 of 13 

36. (New) The computer program product of claim 31, wherein instructions to use 
the structural information to identify a binding region of the protein comprises instructions to: 

map the empty volumes available for ligand binding in the protein to identify one or more 
potential binding regions; 

generate initial conformations for one or more ligands known to bind the protein using 
docking techniques in each of the one or more potential binding regions and scoring a 
preliminary energy function for at least some of the intial conformations; 

select from the initial conformations for each of the known ligands a plurality of best 
conformations in each of the potential binding regions based at least in part on the preliminary 
energy scores; and 

identify the probable binding site based on a spatial location of the best conformations. 

37. (New) The computer program product of claim 31, wherein: 
the annealing molecular dynamics includes a full atom force field. 

38. (New) The computer program product of claim 31, wherein: 
the solvation effects include a continuum description of solvation. 

39. (New) The computer program product of claim 31, wherein: 
the solvation effects include a surface-area based solvation model. 

40. (New) The computer program product of claim 31, wherein 

instructions to calculate a binding energy for ligand in the set of ligands includes taking 
the difference in the ligand energy in the protein and in solution. 

41 . (New) The computer program product of claim 3 1 , wherein: 

the binding energy is calculated for a ligand according to a scoring function comprising 
substracting the free energy of the ligand in water from the energy of the ligand in the protein. 



Applicant ; Wely B. Floriah , Nagarajan Vaidehi, Attorney's Docket No.: 0661S-607002 /CIT3I92-C 

William A. Goddard, 01 

Serial N . : 10/010,725 

Filed : November 30, 2001 

Pag : 12 of 13 

42. (New) The computer program product of claim 31, wherein: 

the binding energy is calculated for a ligand according to a scoring function comprising 
subtracting the free energy of the protein and the free energy of the ligand from the free energy 
of the ligand in the protein. 

43 . (New) The computer program product of claim 3 1 , further comprising 
instructions to: 

identify from the set of ligands one or more iigands predicted to have high binding 
affinity based on the calculated binding energy of the ligands in the binding site. 

44. (New) The computer program product of claim 3 1 , wherein: 

generating and ranking initial conformations includes scoring an energy function for at 
least some of the initial conformations and ranking the initial conformations based at least in part 
on the energy scores. 

45 . (New) The computer program product of claim 3 1 , wherein: 
generating and ranking initial conformations includes determining a percentage of 

the ligand surface area buried within the protein for each of the initial conformations and 
determining energy scores only for a subset of the preferred conformations, wherein each of the 
preferred conformations in the subset has a calculated percentage of the ligand surface area 
buried within the protein which exceeds a predetermined threshold. 



